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ABSTRACT 

We consider the effect of compact stellar remnants on the interstellar medium of 
a massive star cluster following the initial burst of star formation. We argue that 
accretion onto stellar-mass black holes is an effective mechanism for rapid gas depletion 
in clusters of all masses, as long as they contain progenitor stars more massive than 
> 50 M,q. This scenario appears especially attractive for the progenitor systems of 
present-day massive globular clusters which likely had masses above M > 10 7 Mq. 
In such clusters, alternative mechanisms such as supernovae and stellar winds cannot 
provide a plausible explanation for the sudden removal of the primordial gas reservoir 
that is required to explain their complex chemical enrichment history. 

In order to consider different regimes in the rate of gas accretion onto stellar mass 
black holes, we consider both the Bondi-Hoyle approximation as well as Eddington- 
limited accretion. For either model, our results show that the cluster gas can be sig- 
nificantly depleted within only a few tens of Myrs. In addition, this process will affect 
the distribution of black hole masses and, by extension, may accelerate the dynami- 
cal decoupling of the black hole population and, ultimately, their dynamical ejection. 
Moreover, the timescales for gas depletion are sufficiently short that the accreting black 
holes could significantly affect the chemistry of subsequent star formation episodes. 

The gas depletion times and final mass in black holes are not only sensitive to 
the assumed model for the accretion rate, but also to the initial mass of the most 
massive black hole which, in turn, is determined by the upper mass cut-off of the 
stellar initial mass function. Given that the mass function of "dark" remnants is a 
crucial parameter for their dynamical ejection, our results imply that their accretion 
history can have an important bearing on the observed present-day cluster mass-to- 
light ratio. In particular, we show that the expected increase of the upper mass cut-off 
with decreasing metallicity could contribute to the observed anti-correlation between 
the mass-to-light ratio and the metallicity of globular clusters. 

Key words: globular clusters: general - accretion, accretion disks - black hole physics 
- stellar dynamics - stars: abundances - stars: formation. 



1 INTRODUCTION 

Over the last few years, observational evidence has accu- 
mulated that few, if any, globular clusters (GCs) consist 
of a single stellar population with a well-defined age and 
chemical composition. Instead, multiple stellar populations 
have been observed in a number of mas sive Galactic GCs 
(e.g. iGratton. Carretta fe Bra gaglia 201^). Photometric ev- 
idence for multiple populations comes in the form of split 
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sequences in the colour-magnitude diagram (CMD) during 
at least one of the various stellar evolutionary phases, i.e. 
mai n sequence, sub-gi ant branch, horizontal branch, etc. 
(e.g. IPiotto et al]|2007f ). More detailed spectroscopic studies 
have shown that the different stellar populations within the 
same GC often show differences in the relative abundances 
of light elements involved in the proton- capture proce s s such 
as C N, O, F, Na. Mg. Al. or Si (e.g. [5sboriilll97ll : iKraftj 
ll979tlGratton. Sneden fc Carrettall2004 l. 

This can only be explained if the different stellar pop- 
ulations are born from gas with different chemical composi- 
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tions. A popular scenario to explain the time- variable chem- 
istry of the gas reservoir invokes the "pollution" of the pris- 
tine gas by the chemically enriched winds from asymptotic 
giant branch jAGB ) stars of the first stellar generation (e.g 



Gratton et"al]|200ll ; iRamirez fc Cohen|[20o3 ; ICarretta et all 
2010). 



This "self-enrichment" scenario, however, is complicated 
by the fact that the spectroscopic data also demonstrate 
that the different stellar populations in a GC often show 
the same abundances of iron peak and other elements that 
are produced in type II supernovae (SNe II). The chemi- 
cal peculiarities also cannot be limited to surface contam- 
ination, or the y would be diluted after the first d redge- 
up episode (e.g. iGratton. Carretta fc Bragag lia 2012). Cur- 
rently, no scenario offered to explain the presence of mul- 
tiple populations in GCs is able to fully explain the pe- 
culiar abundance anomalies, at least not without resort- 
ing to a non-canonical initial mass funct ion (e.g. iBekkil 
l201ll ; IGratton, Carretta fc Bragaglial [2012). For instance, 
the merging of primordial clusters or gas clumps struggles to 
simultaneously explain both the similarities in the a-element 
abundances and the diff erences in the light element abun- 
dances (e.g. Bekki 2012). The simplest solution to explain 
these peculiar abundance patterns is that an efficient mech- 
anism for gas depletion must have operated after the forma- 
tion of the first generation, but before th e birth of the s econd 
generation (e.g. lConrov fc SpergejlioTll; IConrovll2012l ). 

Building on previous work (ICottrell fc Da Costalll98ll ; 
ISmithlll987l ; ICarretta et al.ll2010l ). the following picture for 
the formation of multiple stellar populations in mass i ve GC s 
has therefore been proposed by IConrov fc Spergell l|201ll ): 
the most massive members of the first stellar generation 
quickly reach the end of their lifetime and explode as SNe 
II. The energy injected into the ISM from these explosions 
can remove the material that is enriched in a-process ele- 
ments from the GC, and can interrupt the star formation. 
Over time, a second gas re servoir is then created from the 
ejecta of either A GB stars (iDecressin et al.ll2007l ) or rotat- 
ing massive stars l|Ventura et al. 200 ll ). combined with any 
remaining (and/or additional) pristine gas. After a few hun- 
dred Myrs, a second generation of stars is born from this 
mix of processed and unprocessed material. To first order, 
this scenario could explain the observed abundance varia- 
tions in the light elements while avoiding differences in the 
a-element abundances. 

However, the question is whether the SNe II can in- 
deed remove (only) the gas component that is enriched 
in a-elements. Theoretical arguments suggest that this 
is only feasible in clu s ters less massive than ~ 10 7 M 
jDopita fc Smith! Il98fj ; iKrause et all l2012h . because more 
massive clusters contain enough gas mass that the SNe 
II cannot plausibly unbind it all. The masses of present- 
day GCs are likely to be only a small fraction of their 
original mass (e.g. iD'Ercole et alj 120081 ) . It follows that 
for the most massiv e Galactic G Cs, including u j Cen (e.g. 
iPancino et afll2003h. NG C 1851 jBekki fc Yond 120121) . 4 7 
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2012al) . NGC 6752 (|Milone et al.l 120101. 



M22 1 Milone et al.l2012bl ). and others, it appears likely that 
the initial gas reservoir had a mass of at least ~ 10 7 Mq, and 
probably more. It thus appears doubtful that SNe II were 
able to efficiently remove the a-enhanced gas from these 
clusters. 



A number of alternative scenarios for the "purging" of 
the gas reservoir in massive stell ar clusters ca n be f ound 
in the literature. For example, iHerwig et al. (2012) in- 
voke Galactic plane passages and the associated ram pres- 
sure stripping of the gas content to explain the enrich- 
me nt history of ui Cen, t hus elaborating on earlier work 
bv lTavler fc Wool |l97l). 

However, in order for this sce- 
nario to explain the multiple stellar populations in an ar- 
bitrary GC, some fine-tuning between their age difference 
and the orbital period of the cluster may be required. In 
many cases, the orbital period could be altogether too long 
for this mechanism to be viable, particularly for clusters 
at large Galacto-centric radii. Another alternative was pro- 
posed by ISpergell |l99ll) who invokes the relativistic winds 
of millisecond pulsars as a catalyst for gas expulsion. Lastly, 
the feedback fro m x-ray binary jets onto the ISM cou l d also 
contribute (e.g . Fender, M accarone fc van Kesterenl 120051 ; 



lJustham fc Schawinskill2012l ) 

Although the precise mechanism for gas depletion in 
primordial GCs is currently unknown, the effects on the 
sub sequent cluster evolution can be dr amatic. For exam- 
ple, [Marks^^rouEa^&^aumgard^ (2008) showed that rapid 
gas expulsion can unbind a significant fraction of stars in 
the cluster outskirts. It can even result in complete clus- 
ter dissolution, e specially in the absence of primordial mass 
segregation (e.g. iMarks fc Kroupal 2010). Thus, in addition 
to accounting for the peculiar abundance anomalies, gas de- 
pletion could play a crucial role in deciding the present- 
day relative numbers of first and seco nd generation stars in 
GCs hosting multiple populations (e.g. ID'Ercole et al.ll2008l ; 
IConrovllioia ). 

In this paper, we consider yet another mechanism for 
the removal of gas after an initial burst of star formation, 
namely the accretion of the ISM onto compact stellar rem- 
nants such as black holes and neutron stars. This scenario is 
appealing because it can significantly deplete the gas reser- 
voir even for very large initial cluster masses, and does not 
require any specific orbital conditions. We will show that, 
in principle, the most massive stellar black holes with par- 
ent stars at the high end of the initial mass function (IMF) 
can be very effective in depleting the gas within the cluster. 
This can significantly modify the final black hole (BH) mass 
distribution, and therefore accelerate the phase of BH-BH 
ejections that, as theoretical st udies suggest, should remove 
most BHs from the c luster (e.g. lPhinnev fc S igurdssor] |l99ll ; 
iDowning et al.ll2010l ). 

The high-mass end of the stellar initial mass function 
of primordial GCs and, by extensi on, the numb er of BH 
progenitors is highly uncertain (e.g. M ae dcr 2009])- In open 
clusters and star-forming associations, however, it has been 
shown that the maximum IMF mass correlates with the total 
cluster mass, and stars sufficiently massive to be the progen- 
itors of BHs should form in clusters more massive than ~ 10 3 
M (e.g. lKirk fc Mversll201ll . l2012l ). Thus, we expect that at 
least some BHs should have formed in a typical primordial 
GC. In support of this, there now exists compelling observa- 
tional evidence in favou r of BHs residing in present-day GCs. 
iMaccarone et all (|2007l ) confirmed the first such BH in the 
giant elli ptical galaxy NGC 4 472 in the Virgo Cluster. More 
recentlv. lStrader et al.l (|2012l ) reported two flat-spectrum ra- 
dio sources in the Galactic globular cluster M22, which the 
authors argue are accreting stellar-mass BHs. This result 
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suggests that the dynamical ejection of BHs from clusters 
may not be as efficient as previously predicted, and that a 
cluster such as M22 could contain on the order of ~ 5 — 100 
stellar-mass BHs. Observational evidence has also been re- 
ported in favour of an intermediate- mass black hole (IM BH) 
existing in the globular cluster M54. Ilbata et al.l (|2009T l first 
reported the detection of stellar density and kinematic cusps 
in M54, which they interpreted as being due to the pres- 
ence of a central BH with a mass of ~ 940 Mq . More 
recently, however, IWrobel, Greene fc Hoi (|201lh used Chan- 
dra and Hubble Space Telescope astrometry to rule out the 
existence of an x-ray counterpart to the proposed IMBH, 
and placed an improved upper limit on its luminosity of 
L(8.5 GHz) < 3.6 x 10 29 erg s" 1 . Nevertheless, the observa- 
tional evidence supports the presence of at least some BHs in 
present-day GCs. Assuming that BHs are efficiently ejected 
during the dynamical evolution of the cluster, it seems likely 
that a more substantial population of stellar-mass BHs once 
existed in primordial GCs. 

Following this introduction, we describe in §[2] our as- 
sumptions for the initial conditions, i.e. the stellar mass 
function and dynamical state of the cluster at the time of BH 
formation. Using an analytic closed-box model for the evo- 
lution of the remnant mass function due to accretion from 
the ISM, we consider two theoretical models for accretion 
of ISM onto a black hole, namely the Bondi-Hoyle approx- 
imation and the Eddington limit. These models, which are 
described in more detail in §[3] represent two different mass- 
dependencies for the accretion rate. The results of our analy- 
sis, which are presented in §|4] show that the accretion of gas 
onto compact stellar remnants can have a significant effect 
on subsequent episodes of star formation, and even present- 
day mass-to-light ratios. We discuss the implications of our 
findings on models of GC formation in §0 and summarize 
our work in §[6j 



We assume an initial mass function for the stellar pop- 
ulation according to: 



t l \ dN R - 



(1) 



where N is the number of stars with a given stellar mass m, 
and a a nd (3 are consta nts. We assume a power-law slope of 
a = 2.3 l|Salpetedll955h . This choice minimizes the number 
of high-mass stars relative to other IMFs used in the litera- 
ture, so that the estimates we derive in §|4] for the number 
of stellar remnants and the amount of gas accreted by them 
are conservative. 

As a technical point, the parameter j3 is needed to en- 
sure that the correct total stellar mass is preserved when 
integrating the IMF. It is determined by normalizing Equa- 
tion [1] using the total stellar mass: 



AL 



f m (rn)mdm, 



(2) 



where m m in and m m ax are the minimum and maximum stel- 
lar masses, respectively. Integrating Equation[2]with respect 
to m and solving for j3 yields: 



M s (2-a) 



(3) 



Plausible values for the upper and lower mass cut-off at 
low metallicity are m m j n = 0.0 8 Mq and m max = 150 Mq 
fe.g. iDabringhausen et al]|201ll ). However, we consider sev- 
eral different values for m max in order to test the sensitivity 
of our results to this assumption. The mass cut-off values 
are chosen to be conservative for very low metallicity since, 
as we will show, a higher initial maximum BH mass trans- 
lates into a shorter gas depletion timescale. It follows that 
our results correspond to upper limits for the gas depletion 
timescale due to accretion onto stellar-mass BHs. 



2 INITIAL CONDITIONS 

In this section, we present an analytic model for the forma- 
tion of a massive star cluster from its parent giant molecular 
cloud. 



2.1 Progenitor Mass Function 

We consider a molecular cloud of total mass M that con- 
denses to form a star cluster at time t = 0. We assume 
a star formation efficiency (SFE) of 25%, so that, after 
the cluster is born, the total mass in stars and gas are 
M s = 0.25M and M g = 0.75Af, respectively. While this 
value for the SFE agrees with availa ble theoretical con- 
straints (e.g. iMcKee fe Ostrikerl l2007h . we point out that 
our results are insensitive to the precise value of the SFE, 
as long as it is neither extremely low or extremely high. In 
the first limit, the stellar mass of the infant cluster would 
be too small for many massive stars to form, and hence the 
number of BHs would be limited. At the other extreme, an 
extremely high SFE would leave little gas, and hence the 
BHs could not grow significantly in mass. Our chosen value 
for the SFE stays safely away from these regimes. 



2.2 Progenitor Velocity Distribution 

In order to calculate the relative velocity between the accre- 
tor and the gas (which affects the accretion rate) as well as 
the retention fraction of stellar remnants after they experi- 
ence a potential natal kick, the velocity distribution of the 
progenitor population is required. 

In this respect, we assume that the core radius r c is 
initially comparable to the half- mass radius i\, since pri- 
mordial clusters are expected to be more extended than 
their (dynamically evolved) prese nt-day counter-parts (e.g. 
ISpitzejfl987l : iHeggie fc HutH2003T ). We further assume pri- 
mordial mass segregation, i.e. that all remnants form inside 
r c , and that energy equipartition has been achieved within r c 
(we will come back to these assumptions in §[5|. We therefore 
adopt a Maxwell-Boltzmann velocity distribution for the 
proge nitors at every stellar mass (e.g. iBinnev fc Tremainel 
1 19871 ): 

f v (m,v) = N(m)e-? v2/ ° (m)2 , (4) 

where N(m) is the number of stars with mass m, <r(m) is the 
velocity dispersion for stars of mass m, and v is the stellar 
velocity. For the velocity dispersion <j(m), we assume 
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where ao is the central velocity dispersio n of the average stel- 
lar m ass m. For a King model, this is l|Binnev fe Tremaind 
Il987t) : 



UirG,_ 
(To = (^-g-(Ps + I 



1/2 



(6) 



where r c is the core radius, p s is the average stellar mass 
density in the core, and p s is the average gas mass density 
in the core. 



2.3 Progenitor Lifetimes 

For the progenitor lifetime r p , we assume that the main- 
sequence (MS) lifetime TMsijn) provides a good approxima- 
tion, provided the progenitor mass is m 18 M©. This 
seems justified because the MS lifetimes of low-mass stars 
greatly exceed that of every other evolution ary phase, typ- 
ically by several orders of magnitude (e.g. IClavtonl 1 19681 ; 
Ilbenlll99ll ; lMaeder| [2009) . Note that TMs(f7i) is a function of 
only the stellar mass. In other words, we are ignoring any 
metallicity-dependence, since metallicity should only weakly 
affect the total stellar lifetime. For the MS lifetime, we as- 
sume (|Hansen fc Kawaler]|l994h : 



TM3 (m) = T (^-) 



(7) 



with to = 10 10 years. For progenitor masses m > 18 Mq , we 
impose a fixed total lifetime of 7 Myrs (note that Equation]?] 
yields the same MS lifetime of 7 Myrs for m = 18 M Q ). 
This is in rough agreement with stellar evolution models 
(Rob Izzard, private communication), which predict a near 
constant l i fetime for massive stars at low metallic ty (e.g. 
llbenlll99ll : lHurlev. Pols fc Toutll200ol ; lMaedej|200gh . Thus, 
our final estimate for the total progenitor lifetime is: 



TVi = max 



(T A fs(m),7Myrs). 



(8) 



2.4 Initial Remnant Mass Function 



For th e initial remnant m asses, we use Equations 4, 5, 7, and 
9 from I Fryer et al.l ((2012) for progenitors in the mass ranges 
7.2 < m/M < 11, 11 < m/M < 30, 30 < m/M Q < 50, 
50 m/M Q < 90, and m > 90 Mq, respectively. In that 
order, the initial remnant masses are: 



M r 



1.35M, 



0i 



(9) 



1.1 + (Ue^- 11 ' '/ 4 ' - 2.01e°' 4(M p- 26 ' 0) , (10) 



M r = min(33.35+4.76(Af p -34.0), Af p -0.1(1.3M P - 18.35)), 

(11) 



and 



Mr = 1.8 + 0.04 x (90.0 - M p ), 



M T = 1.8 + log 10 (M p - 89.0), 



(12) 



(13) 



where M p and M T are the progenitor and remnant masses, 
respectively, in units of solar masses, and we have used 
Z = O.Ol^o, with Zq denoting solar metallicity. This choice 
of metallicity roughly agreees with the average for Milky 
Way GCs (e.g. lHarris!ll996. 2010 updated , and should be ap- 
propriate for a star cluster composed of Population II stars. 



For each progen itor mass range, two equations are given 
in lFrver et al.l |2012i ') for the initial-final mass relations, one 
assuming a delayed supernova explosion and the other as- 
suming a rapid explosion. The delayed explosion model pre- 
dicts considerably greater kick velocities that are roughly 
consis tent with what is o bserved for pulsars in the Galactic 
field (|Hobbs et al.ll2005h . Therefore, we assume a delayed 
supernova explosion in choosing our initial-final mass rela- 
tions. 

However, we note that our results are not sensitive to 
this choice. This is because the gas accretion is completely 
dominated by the most massive BHs, specifically those with 
masses > 50 Mq. In this high- mass regime, the differences 
in the final remnant mass functions be t ween the delayed and 
rapid explosion models of iFrver et al.l (|2012h are negligible. 
For the same reason, we do not discuss gas accretion onto 
NSs and white dwarfs (WDs) further: they do not affect the 
results presented here. Nevertheless, in order to check the 
plausibility of our assumptions, we compare the retention 
fractions of NSs to previous estimates in the literature in 

§m 

2.5 Remnant Velocity Distribution 

The accretion rate depends on the relative velocity between 
the gas and the accretor. Thus, it is important to obtain a 
realistic final remnant velocity distribution from the initial 
progenitor velocity distribution. 

The primary factor that determines the initial remnant 
velocity distribution is natal kicks. These are imparted to 
most NSs, and some (low-mass) BHs, at the time of their 
formation. The precise physical mechanism responsible is 
unknown, but appears to be related to asymmetries in the 
collapse of the progenitor core, or the subsequent supernov a 
explosion (e.g. IPfahl. Rappaport fc Podsiadlowskil l2002al ). 
Typical kick speeds range from 100 to 200 km s _s , although 
radio pulsars wi th speeds in excess of 10 km s -1 have been 
observed (e.g. Lvne fc Lorimerl 1 19941 : ICordes fc Chernofil 
1 19981 : iHobbs et al.1 120051 1. and kick speeds < 50 k m s" 4 
are also t ho ught to be possible (e.g. iBlaauwl Il96l | : 
Pfahl et all l2002bl: [Pfahl. Rappaport fc Podsiadlowskil 



2002bl; iPodsiadlowski et all 120041: ISmits et all \20m. 



Paolillo et al.ll201lD . 



We adopt different assumptions for the kick speed 
for different progenitor mass ranges. These assumptions 
are based on the physical processes thought to be driving 
the kicks. Here, we present only our assumed kick veloc- 
iti es for each mass r ange, and refer the interested reader 
to IFrver et all (|2012l ) for more details. Specifically, we as- 
sum e that no kick occurs for progenitor masses < 11 Mq 
(e.g. [Podsiadlowski et al. 2004! ). For progenitors in the mass 
range 11 < m/M < 30 (|Frver fc Kalogera|[200lT ). we adopt 
a kick velocity Vkick by randomly samp ling from the distri- 
bution of observed velocities reported bv lHobbs et al.l (|2005l ) 
for radio pulsars in the field of our Galaxy. For progenitors 
with masses 30 ^ m/M© < 40, we adopt a constant kick ve- 
locity of 50 km s~ , since the amount of fallback is thought 
to be less sensitive to the explosion mechan ism in this range 
of progenitor masses (e.g. lFrver et aljboi/t ). Finally, for pro- 
genitors with m/M© ^ 40, we assume prompt collapse to a 
BH without a supernova explosion, and thus without any 
kick. 
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In order to calculate the initial velocity distribution for 
our remnant population, we add the kick velocity Ukick to 
a randomly drawn velocity vo from Equation [3] for the ap- 
propriate progenitor mass. This is done using randomized 
vector addition, since there is no reason to expect that the 
direction of the kick should be correlated with the direction 
of motion of the progenitor through its host cluster. The re- 
sulting velocity w re j is then compared to the escape velocities 
from both the core and the cluster. The latter quantities are 
calculated using the relations 



and 



2GM C 



2GM 



(14) 



(15) 



where M c is the total (stellar and gas) mass of the core, M 
is the total (stellar and gas) cluster mass, and r c is the core 
radius (recall that we have assumed that r c m r& initially). 

A comparison between v Te \ and v esc determines the NS 
and BH retention fractions immediately after their forma- 
tion. If the final velocity is less than the escape velocity of the 
core (i.e. v IC \ < « csc ,c), we leave it unchanged from the ve- 
locity previously drawn from the distribution in Equation [4] 
This is because we assume that the remnant will rapidly re- 
achieve energy equipartition with the rest of the core, where 
the local relaxation time is much shorte r than the global 
relaxation time (e.g. iHeggie fc Hutl [20031 ). If the remnant 
ends up with a velocity greater than the escape velocity 
from the cluster immediately after it receives its kick (i.e. 
Wrel > feac), then it is removed completely from the remnant 
population. However, if the post-kick velocity exceeds the 
escape velocity from the core but not that from the cluster 
(i.e. i>csc,c < Vrel < Ugsc), then the remnant will temporarily 
leave the core but still remain bound to the cluster. In this 
case, we calculate the radius to which the remnant is kicked 
(called the "kick radius") using conservation of energy. This 
requires two additional assumptions, namely that the kick 
is always imparted at a distance r c /2 from the cluster cen- 
tre, and the resulting orbit of the kicked remnant within its 
host cluster is entirely radial. We then calculate the fall-back 
time tfb, which decides when the remnant will return to the 
core. This is given by: 



tfb 



2v t , 



GMc 



(16) 



3 ACCRETION MODELS 

The Bondi-Hoyle approximation describes spherically sym- 
metric accretion, and is based on the assumption that the 
forces due to gas pressure are insignific ant compared to grav- 
itational forces |Bondi fc Hovlelll944h . The background gas 
is treated as uniform and either stationary or moving with 
constant velocity relative to the accretor. This assumption 
gives reasonable accretion rates, provided the properties of 
the gas are such that the density and total angular mo- 
mentum are low, at least in the vicinity of the accretor. 
Numerical studies have confirmed that this is indeed the 
case, provided t he gas is not moving rapidly relative to 
the accretor (e.g. iFrvxell fc Taamlll98cl : lRuffert|[l99i . Il997l : 
iFoglizzo fc Ruffert 19991 ). High accretion rates at least ap- 
proaching the Bondi-Hoyle rate hav e also been observed in 
several high-mass x-ray binaries (e. g. lAgol fc Kamionkowskil 
l2002l : lBarnard, Clark fc Kolbll2008h . It thus appears that, at 
least in the low-density, low-angular momentum regime, the 
Bondi-Hoyle approximation is not unrealistic. 

The Bondi-Hoyle approximation should be regarded as 
a strict upper limit to the true accretion rat^E F° r large 
accretor masses, the Bondi-Hoyle rate specified in S I3.1l can 
become extremely high, and pressure forces could play an 
important role in reducing the accretion rate. Notwithstand- 
ing, theoretical studies have confirmed that super-Eddington 
accretion is indeed possible for massive BHs. For example, 
although the accretion r ates of SMBHs at high reds hift are 
poorly constrained (e.g. iHopkins fc Quataerq l2010t ). some 
have bee n observed at redshift z ~ 6 with masses > 5 x 10 9 
M (e.g. lBarth et al.ll2003l : lwillott, McLure fc Jarvisll2003h . 
The prese nce of such massive BH s only « 10 9 years after the 
Big Bang (|King fc Pringlell2006l ) is difficult to explain with- 
out invoking accretion rates approaching the Bondi-Hoyle 
limit. This should be kept in mind when discussing accret- 
ing BHs in primordial globular clusters, which were also 
born soon after the Big Bang and at similarly low metal- 
licity. Additionally, photon-trapping can occur at very high 
accretion rates, which makes accretion disks radiatively inef- 
ficient and_pj^y2des_ain^ns_£fcircumventing the Eddington 
limit jPaczvnskv fc Wiitalll980l ). Otherwise, when an accre- 
tor is radiating at above the Eddington luminosity, signifi- 
cant amounts of gas can be expelled at hig h velocities due 
to th e intense winds that are initiated (e.g. iKing fc Pounds! 
l2003h . 

Eddington-limited accretion applies when pressure 
forces become important, and the outward continuum 
radiation force balance s the inward gravitational force 
(Eddingtonl ll92ol [l~930h The accretion luminosity corre- 



This assumption is justified given that primordial GCs are 
thought to have been more extended at birth than they are 
now. It follows that stars should be relatively unaffe cted by 
dynamical friction once outside t he core (e.g. IMapelli et al.l 
l2006l : lLeigh. Sills fc Kniggell201 ll ). 

If tms +tfb < Tf , where n denotes the time from cluster 
birth until the gas reservoir is depleted, then the remnant 
will return to the core with enough time to accrete at least 
some gas from the dense ISM. Once in the core, the remnant 
quickly re-achieves energy equipartition with the surround- 
ing population. The remnant's final velocity is then once 
again given by Equation [4] 



1 In this context, it is sometimes useful to distinguish between 
a capture rate and an accretion rate. The former describes the 
rate at which material becomes bound to the accretor, at least 
temporarily, whereas the latter describes the rate at which the 
accretor actually grows in mass. The Bondi-Hoyle approximation 
can provide a good estimate for the capture rate, whereas the true 
accretion rate may well be much lower. For the purpose of this 
paper, however, we assume that all captured gas is accreted by 
the BH, which is another argument for regarding it as an upper 
limit. 
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sponding to this limit is (|Rvbicki fc Lightmarj|l979h : 

4nGmc 

t-iKdd = . 



(17) 



where m is the accretor mass, c is th e speed of light, and k i s 
the electron scattering opacity fe.g. iKing fc Pounds! [20031 . 
Compared to the Bondi-Hoyle approximation, the Edding- 
ton limit provides a much more conservative estimate for the 
accretion rate in the limit of high gas density. This is due 
the linear dependence on the accretor mass m (see § 13. ip . 
Additionally, the Eddington rate should be much closer to 
the true accretion rate if the gas contains significant angular 
momentum, and accretion proceeds mainly via angular mo- 
ment um re-distribution within a disk (e.g. IKing fc Pounds! 
l2003h . 

We assume a constant time-independent background 
density and velocity for the accreting gas, since a more so- 
phisticated treatment is beyond the scope of this paper. To 
illustrate the importance of accreting BHs for prolonged star 
formation in primordial GCs, we use the Bondi-Hoyle and 
Eddington approximations to explore two different regimes 
in the gas properties. The Bondi-Hoyle approximation is rep- 
resentative of the gravity-dominated, low-angular momen- 
tum regime, whereas the Eddington approximation is rep- 
resentative of the pressure-dominated or high-angular mo- 
mentum regime. 

3.1 Calculating the Accretion Rates 

The Bondi-Hoyle accretion rate is given by Equation 2 of 
M accarone fc Zurekl (|2012|). which has been re-scaled from 
the formula of Ho. Terashima fc Okaiimal (|2003[ ): 



Mb 



7 x 10" a M 



©yr 



Vci + v 2 



\Mq) \10 6 cm- 3 / Vl0 6 cms 
(—)' 



A 



(18) 

where m is the accretor mass, n is the central gas density, 
c s is the sound speed, and v is the relative velocity between 
the gas and accretor. We assume n ~ 10 6 cm -3 with central 
temperatures of a few thousand K, which translates into 
c s ~ 10 6 c m s" 1 . These values were arg ued bv lD'Ercole et all 



l|2008h and|M 

accarone fc Zurekl ([2012) to be representative 
of the cores of primordial GCs, yielding A w 10~ 9 Mq yr -1 
for small accretor masses and velocities. 

Th e Eddington-limited a ccretion rate, in turn, is given 
by (e.g. lKing fc Pounds! [2001 ): 

4nGm ,_„. 
M Bd d = , 19 

T]KC 

where c is the speed of light, 77c 2 is the accretion yield from 
unit mass, and k, is the electron scattering opacity. We adopt 
T7 = 0.1 for the accretion efficiency. 



3.2 Evolving the Remnant Mass Function 

The most massive stars in the cluster are the first to reach 
the end of their lives due to stellar evolution. If gas is still 
present in significant quantities, the most massive stellar 
remnants (i.e. BHs) should begin accreting from the sur- 
rounding ISM almost immediately after their formation, 
since they do not experience kicks. 



If a kick is imparted, there is a possibly indefinite delay 
before significant accretion begins, because the remnant can 
be expelled from the core. However, if the remnant remains 
bound to the cluster and the fall-back time tfb is sufficiently 
short, then the remnant will eventually begin accreting from 
the ISM after it returns to the central cluster regions. Thus, 
accretion is assumed to occur only in the core, and to be- 
gin after a total time n — tms + tfb , where tfb = if the 
remnant never leaves the core. Once returned to the core, 
we assume that the remnant remains there. This is not un- 
reasonable since the timescale for accretion to occur is com- 
parable to the crossing time of the core, even for accretion 
rates lower than the Eddington rate. Therefore, the accretor 
velocity should be significantly reduced within a single core 
crossing time. 

The post-accretion mass of a given remnant is calcu- 
lated as follows. The total time spent accreting from the 
ISM can be written as: 



dt. 



(20) 



where the upper limit of integration tj corresponds to the 
time at which the post-accretion remnant mass is calculated, 
and Ti = tms + tfb ■ Re- writing this in terms of the accretion 
rate m = dm/dt gives: 



J m; 



dm 
m ' 



(21) 



where mi and mt are the initial and final remnant masses, 
respectively. We then integrate and solve for mt in order to 
obtain the final remnant mass. Using Equation [18] for the 
accretion rate as an example, this yields the final remnant 
mass: 



mt 



rrii 



1 — m\Ar^. 



(22) 



In our analysis, we have neglected the evolution of 
the remnant velocity distribution during the accretion pro- 
cess even though, in principle, conservation of momentum 
causes the BHs to slow down as they accrete mass. In the 
Eddington-limited case, this does not affect the accretion 
rate, since Eq. [19] is independent of the accretor velocity. 
In the Bondi-Hoyle case, however, a reduction in the accre- 
tor velocity increases the accretion rate, and hence the rate 
of gas depletion. Neglecting the BH deceleration, therefore, 
leads to a lower limit for the accretion rate in the Bondi- 
Hoyle approximation, and properly accounting for these ef- 
fects would only strengthen our conclusions on the efficiency 
of the gas depletion. 

A possibly more important consequence of a lower av- 
erage BH velocity is that it could lead to a more efficient 
migration of massive BHs towards the cluster centre, and 
thus could acc elerate the onset of the dynamical BH ejec- 
tion p hase (e.g. iPhinnev fc Sigurdssor3ll99ll : iDowning et al.l 
2010) or the rate of BH-BH mergers near the cluster cen- 
tre, as discussed in more detail in S I5.1I Both effects caused 
by gas accretion - deceleration and mass growth - should 
therefore be taken into account when constructing dynami- 
cal models of the BH population in gas-rich environments. 
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4 RESULTS 

4.1 Black Hole Depletion Times 

The results of our analysis are summarized in Figures [T] and 
[2]which show the time evolution of the relative mass fraction 
of gas, stars, and remnants for Bondi-Hoyle and Eddington- 
limited accretion, respectively. In both cases, the first BHs 
form 7 Myrs after the birth of their progenitor population, 
and immediately begin accreting from the ISM. 

In the case of Bondi-Hoyle accretion (Figure [T}, the 
steep dependence of the accretion rate on black hole mass 
(M oc m 2 ) causes the most massive BH(s) to reach very high 
accretion rates after only a few Myrs. For example, for an 
IMF upper mass cut-off of 150 Mq (dashed line), the most 
massive BH reaches an accretion rate of M w 10 -4 Mq 
yr -1 after only ~ 4 Myrs. At this point, it enters a phase 
of runaway growth and, as shown in Figure [T] the remnant 
mass fraction rises asymptotically from only a few percent 
to nearly 80% in as little as a few thousand years. Thus, 
in the Bondi-Hoyle case, the entire gas reservoir is depleted 
within a few Myrs of cluster formation, nearly independent 
of our assumption for the IMF upper mass cut-off. It follows 
that the gas depletion time will be under 1 Gyr provided 
the accretion rate > 10 -2 Mb-h ■ 

In the case of Eddington- limited accretion (Figure [2]), 
in contrast, a runaway accretion phase is never reached, be- 
cause of the weaker (linear) dependence of accretion rate 
on BH mass. Consequently, the remnant mass fraction rises 
gradually at first over a period of several tens of Myrs, reach- 
ing 5-10% (depending on the assumed IMF upper mass cut- 
off) within ~ 50 Myrs of cluster formation. Nevertheless, 
the accretion rate continues to rise as the BHs grow in mass 
and, after ~ 100 Myrs, the remnant mass fraction reaches 
20-40%. Thus, provided the accretion rate > 10 _1 MsdcL, ac- 
cretion onto stellar-mass BHs will significantly contribute to 
depleting the available gas reservoir within a few hundred 
Myrs. 

The key conclusion from these plots is that, in both 
accretion models, stellar-mass BHs cause a non-negligible 
reduction in the total mass of the gas reservoir within as 
little as a few tens of Myrs. This is the case even if we 
truncate the IMF upper mass cut-off at 60 Mq, which is 
considerably lower than what is theoret ically pred i cted for 
mass ive stars at very low metallicity (e.g. Iben 1991; Maederl 
12009(1 . 



Our results show that the assumed model for the accre- 
tion rate and its dependence on accretor mass can have a 
significant effect on the final distribution of remnant masses. 
Interestingly, in the Bondi-Hoyle case, the accretion rates 
can become sufficiently high that, if one or a few BHs reach 
these rates well before their peers, there will be a large dis- 
crepancy between the final masses of these BHs and the 
rest. It is even feasible that one especially massive BH (or 
intermediate-mass BH) will completely dominate the accre- 
tion "race". The possible implications of this case for the 
subsequent dynamical evolution of the BH population and 
the evolving cluster mass-to-light ratio will be discussed in 
Section [5] 

Finally, kicks have a minor effect on BH population size 
due to the fact that it is only low-mass BHs that receive 
kicks. For initial cloud masses M = 10 6 , 10 7 , and 10 8 M , 
we find BH retention fractions of roughly 66%, 90%, and 
96% respectively. Kicks also do not significantly affect the 
gas depletion times since low-mass BHs have the lowest ac- 
cretion rates. This assumes, however, that the IMF upper 
mass cut-off is sufficiently large that the most massive BHs 
will not receive kicks (see §[23} . 

4.3 Neutron Stars 

Due to their relatively small masses, accretion leaves the 
numbers of NSs largely unchanged for up to ~ 100 Myrs. 
Even in the Bondi-Hoyle case, many of the first NSs to form 
will only accrete enough material within this time-frame to 
end up as especially massive NSs, as opposed to low-mass 
BHs. Thus, the primary factor affecting NS population size 
in the core is natal kicks. We can therefore consistently com- 
pare the NS retention fractions we obtain for our model 
to previous estimates given in the literature calculated us- 
ing more sophisticated methods. This is meant to provide 
a test of the validity of our assumptions. For initial cloud 
masses M = 10 6 ,10 7 , and 10 s M Q , we find NS retention 
fractio ns of roughly 1 4%, 3 7%, and 73%, respectively. Re- 
cently, llvanova et al.l ()2008h used Monte Carlo models for 
GC evolution to look at NS retention due to natal kicks. 
Our results are in close agreement with theirs for an initial 
cloud mass M = 10 6 Mq with 25% star formation efficiency, 
which yields a total stellar mass that agrees roug hly with the 
cluster mass considered by llvanova et al.l <|2008h . 



4.2 Black Hole Growth 

The effects of a larger initial BH mass are also demonstrated 
in Figure [3l which shows the maximum BH mass as a func- 
tion of time for three different assumptions for the upper 
mass cut-off for the IMF, namely 150, 100, and 60 M Q . In 
the Bondi-Hoyle case, the upper limit for each mass is ob- 
tained by dividing the remaining gas mass by the number of 
BHs in the highest-mass bin (2, 6, and 19 for an IMF upper- 
mass cut-off of 150, 100, and 60 Mq , respectively) at the time 
when runaway accretion sets in. This represents a theoreti- 
cal upper limit to the maximum final BH mass. We caution 
that, in the Bondi-Hoyle case, such extreme BH masses con- 
flict with curre nt observational upper limits in several mas- 
sive GCs (e.g. iMaccarone fc Servillatl 120081; iMiocchil l 20ld ; 
iKirsten fc Vlemmingsll2012l ; iMcNamara et al.ll2012h . 



5 DISCUSSION 

We now discuss the implications of our results for the dy- 
namical ejection of BHs, subsequent episodes of star forma- 
tion, and the present-day cluster mass-to-light ratio. 



5.1 Dynamical BH Ejection 

Several theoretical studies suggest that BHs which form soon 
after the formation of a massive star cluster should decouple 
dynamically from the rest of the stellar population due to 
their much larger masses. Consequently, if they don't already 
form in the cluster core they w i ll rapidly segregate there (e.g. 
Sigurdsson fc Hernquistlll993l ; iPortegies Zwart fc McMillanl 
20001 ). Inside the core, BH-BH binaries form dynamically 
through three-body scattering interactions, and most BHs 
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5 10 15 20 

Time (in Myrs) 

Figure 1. The relative mass fractions in gas (blue), stars (red), 
and stellar remnants (black) as a function of time (in Myrs), cal- 
culated using the Bondi-Hoyle approximation to estimate the ac- 
cretion rate. A star formation efficiency of 25% and an IMF slope 
of a = 2.3 were adopted. The dashed, dotted, and solid lines cor- 
respond to an IMF upper mass cut-off of 150, 100, and 60 Mq, 
respectively. This plot is insensitive to our assumption for the 
total cloud mass M, provided M > I0 7 Mq. 





Time (in Myrs) 

Figure 2. As Fig.[T] but using the Eddington limit to estimate 
the accretion rate. 



are thought to be eventually ejected from the cluster due to 
strong gravitational interactions with other BH-BH binaries. 
However, it is expected to take at least a few Gyrs for the 
bulk of the primordial BH population to be depleted in this 
way l|Q'Learv et al.l [20061 ^1 . and some studies suggest that 
the required time could even exceed a Hubble time for some 




log Time (in Myrs) 

Figure 3. Mass of the most massive BH as a function of time 
for the case of Bondi-Hoyle (red curves) and Eddington-limited 
(black curves) accretion. The dashed, dotted, and solid lines cor- 
respond to an IMF upper mass cut-off of 150, 100, and 60 Mq, 
respectively. 



clusters dDowning et al.ll2010T i. Recent studies also suggest 
that, in general, dynamical BH ejection may not be effective 
enou gh to remove the entir e BH population fro m the cluster 



(e.g. iMoodv fc Sigurdssonll2009l ; iMackev et al.ll2008T l 



In the absence of BH-BH mergers, we expect accretion 
to accelerate the rate of dynamical BH ejections, because 
mass growth of a BH will increase its gravitationally-focused 
cross-section for encounters. Consequently, if most BHs ex- 
perience a significant mass growth, the average time between 
close encounters (i.e. fly-bys) involving BHs should also de- 
crease. Moreover, the velocity of an accreting BH will be 
reduced significantly if its mass increases on time scales that 
are short compared to the central relaxation time so that it is 
not re-accelerated by the cluster potential. In this case, con- 
servation of momentum predicts that the final post-accretion 
velocity is: 



mi 
1 

m f 



(23) 



where Vi is the initial BH velocity, and mj and mf are the 
initial (pre-accretion) and final (post-accretion) BH masses, 
respectively. Such a reduction in the average BH velocity 
will also reduce the average time between BH-BH encoun- 
ters, which is linearly proportional to the relative velocity 
between interacting objects. For example, an increase in the 
average BH mass by a factor of ~ 10, and the ensuing de- 
crease of the average BH velocity by a factor of ~ 10 could 
result in a decrease in the average time between close BH- 
BH encounters by a factor of ~ 100. This should probably 
be regarded as an upper limit, however, since mergers be- 
tween BHs are likely inevitable at very low BH velocities. 
Nonetheless, we conclude that an episode of accretion from 
the ISM can significantly reduce the BH ejection timescale 
by a factor of several tens or more. 
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5.2 Multiple Episodes of Star Formation 

If accreting BHs are indeed responsible for depleting the 
gas reservoir of a young stellar cluster, a second gener- 
ation of stars can only be born after most BHs are dy- 
namically ejected, because only then can gas re-accumulate 
within the cluster^ Previous studies have shown that the 
amount of gas built up from the winds of evolved stars 
and retained within the cluster c an be significant (e.g. 
iNaiman. Ramirez-Ruiz fcTinl l2012h . and can be sufficient 
for a second generation of stars to be born. Moreover, 
in order to explain the peculiar abundance anomalies ob- 
served in GCs hosting multiple populations, it seems that 
some pristine gas must also somehow re-accumulate within 
the cluster after any polluted material due to SNell from 
the first generation is removed from the gas reservoir (e.g. 
IConrov fc Spergelll201ll ; IConroy||2012i i. 

As discussed in the previous section, our results suggest 
that the epoch of BH ejections stimulated via dynamical in- 
teractions between BH-BH binaries could end much sooner, 
depending on how much the BHs increase their mass via 
accretion. Based on the result s of previous theoretical esti- 
mates for this times cale (e.g. IPhinnev fc Sieurdssonl Il99ll : 
iDowning etafl l201Ch and typical age sprea ds estimated 
for m ultiple populations in massive GCs (e.g. IPiotto et al.l 
we conclude that a scenario in which accreting BHs 
significantly deplete the gas reservoir and are nearly all 
dynamically-ejected from a primordial GC before the birth 
of a second stellar generation is not inconsistent with cur- 
rently available observations and theoretical expectations. 

We caution, however, that the timescales for gas deple- 
tion and BH ejections are somewhat sensitive to the degree 
of initial mass segregation. This is because, at the large ini- 
tial cluster masses thought to be need ed to explain the obser- 
vations (e.g. IConrov fc S pcrgcl 201lj ). the global relaxation 
time (corresponding to the average stellar mass) should be 
comparable to or exceed the age of the Universe. However, 
whenever a remnant forms, its progenitor is the most mas- 
sive star in the cluster, and many BHs should be born with 
masses ~ 100 tim es larger than the average stellar mass (e.g. 
iFrver et alJl2012l ). Therefore, the progenitors of BHs will 
likely migrate into the core first, given that the timescale 
for mass s egregation for a star of mass m is approximately 
l|Vishniadll97Sh : 

iffy 

r scg (m) « — t rh . (24) 
m 

Here, m is the average stellar mass and t T h is the half-mass 
relaxation time (which can be used to estimate the global 
rate of two-body relaxation). Once a star is inside r c , the 
local relaxation time can be shorter than the global or half- 
mass relaxation time by a factor of 100 or more fe.g. lSpitzerl 
Il987t ). The key point is that, with or without primordial 
mass segregation, the stellar remnants of the most massive 
progenitors are likely to be found in the cluster core very 
early on, and to be dynamically relaxed. 

As a mechanism for gas depletion, accretion onto stellar- 
mass BHs has the advantage that it can potentially help to 

2 It is sufficient that they are kicked to sufficiently large clus- 
ter radii that the local dynamical friction timescale exceeds the 
timescale for star formation in the core. 



explain why stars belonging to the second stellar generations 
in massive GCs do not alwa ys show evidence for enrichment 
in iro n-peak elements (e.g. iGratton. Carretta fc Bra gaglia 
2012). This is because the accreting BHs can remove the pol- 
luted material from the SNell of the first generation before 
the second generation is born. However, we have not consid- 
ered dynamical ejections that occur between massive stars 
while still on the main-sequence, an effect that should be 
considered in future studies. Finally, the effects of accret- 
ing BHs on the first episode of star formation should also 
be considered, since it may be important if star formation 
is an extended process. In fact, if the most massive stars 
reach the end of their lives while star formation is still ongo- 
ing, gas accretion onto their BH remnants could reduce the 
number of first generation stars, and affect the shape of the 
IMF. To summarize, independent of the exact scenario, our 
results demonstrate that accretion onto BHs could signifi- 
cantly affect the star formation histories of massive globular 
clusters. 

On the other hand, while the BHs of the first stellar 
generation offer a potentially attractive mechanism for gas 
depletion, their continued gas accretion would pose a prob- 
lem for the formation of the second stellar generation, unless 
the BHs are efficiently removed from the cluster core after 
their formation. This should be considered in any proposed 
scenario explaining the complex star formation history of 
globular clusters. 

5.3 Black Hole Winds 

Depleting the gas reservoir via accretion does not necessar- 
ily require that the accreted matter end up bound to the 
accretor(s). If the accretion rate surpasses the Eddington 
limit, powerful winds could develop that might expel an y 
remaining gas from the cluster (e.g. iKing fc Pounds! [20031. 
It is possible that much of the emission could be in the UV, 
which could in turn result in photoionization heating of the 
ISM beyond the escape velocity from the cluster. Properly 
quantifying the implications of these effects requires detailed 
numerical modeling. This technology is currently unavail- 
able in light of the considerable compu tational difficulties in 
achie ving the necessary resolution (e.g. iHopkins fc Quataerd 
l2010h . 

5.4 The Cluster Mass-to-Light Ratio 

The considerable uncertainties in the accretion rate onto 
BHs, the stellar IMF, and the total amount of accreted ma- 
terial make it difficult to accurately predict the final mass- 
to- light ratio of the stellar cluster from our simplistic model. 

In general, if the BHs accrete from the surrounding ISM, 
and are retained by their host cluster (e.g. the BHs merge 
rather than being dynamically ejected), one may expect a 
correlation between the cluster metallicity and its mass-to- 
light ratio. This is because theoretical studies suggest that 
the IMF u pper mass cut-off increases wi th decreasing metal- 
licity (e.g. (Abel. Bryan fc Normanl2002l ). which, in turn, ac- 
celerates the mass growth of BHs (see §[4]). Interestingly, 
such an anti-correlation between the metallicity of globular 
clusters a nd their mass-to-light rat i o has indeed been ob- 
served bv lStrader. Caldwell fc Sethi l|201ll i in M31. 
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If, on the other hand, higher accretion rates result in 
a more efficient dynamical ejection of BHs, as discussed in 
§ 15.11 such that (nearly) all BHs are dynamically ejected, one 
may even expect a correlation between the cluster metallic- 
ity and mass-to-light ratio. We conclude that a more sophis- 
ticated model is needed to constrain the dependence of the 
final mass-to-light ratio on metallicity. The key point aris- 
ing from our model is that accretion from the ISM onto a 
population of remnants in a primordial GC can significantly 
affect its present-day mass-to-light ratio. 



6 SUMMARY 

In this paper, we have considered accreting black holes as a 
mechanism to deplete primordial GCs of their gas after the 
birth of their progenitor population. This mechanism should 
be effective in clusters of all masses, as long as they are able 
to form massive stars that turn into BHs. 

We have considered both the Bondi-Hoyle approxima- 
tion and the Eddington limit in order to explore two different 
regimes in the mass-dependency of the accretion rate. Our 
results suggest that accreting BHs are able to deplete the 
entire gas reservoir in as little as ~ 10 Myrs, if the Bondi- 
Hoyle rate is assumed. Even at significantly lower accretion 
rates on the order of the Eddington limit, BHs can accrete 
a significant fraction of the total gas reservoir within a few 
tens of Myrs. This timescale is sufficiently short that ac- 
creting BHs would significantly impact the chemistry of any 
subsequent episodes of star formation. The accretion rates 
are sensitive to the IMF upper mass cut-off, in the sense 
that reducing the mass of the most massive BH by a factor 
of ~ 2 will result in gas depletion times that are longer by 
a factor of ~ 2. 

At the end of the accretion period, the average BHs will 
have significantly increased. This is likely to accelerate the 
dynamical decoupling of the BH population evolution from 
the rest of the cluster system, and thus may result in a more 
efficient ejection of the BH population from the cluster. We 
have discussed the implications of this effect on subsequent 
episodes of star formation. 

Lastly, we have pointed out that accretion onto stellar 
remnants may significantly affect the present-day mass-to- 
light ratio, although we are currently unable to precisely 
quantify the magnitude of this effect. 
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